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We estimate the event-by-event (e-by-e) distribution of the ratio (cr) of the magnetic field energy 
to the fluid energy density in the transverse plane of An-An collisions at ^snn = 200 GeV. A Monte- 
Carlo (MC) Glauber model is used to calculate the a in the transverse plane for impact parameter 
b=0, 12 fm at time ti ~0.5 fm. The fluid energy density is obtained by using Gaussian smoothing 
with two different smoothing parameter ag—t).26 , 0.5 fm. For 5 = 0 fm collisions a is found to be 
<C 1 in the central region of the fireball and cr ^ 1 at the periphery. For b=12 fm collisions cr > 1. 
The e-by-e correlation between cr and the fluid energy density (e) is studied. We did not find strong 
correlation between a and s at the centre of the fireball, whereas they are mostly anti-correlated at 
the periphery of the fireball. 


I. INTRODUCTION 

The most strongest known magnetic field (B^ 10^^ — 
10^^Gauss) in the universe is produced in laboratory ex¬ 
periments of An-An or Pb-Pb collisions in the collider 
experiments at Relativistic Heavy Ion Collider (RHIC) 
and at Large Hadron Collider (LHC). Previous theoreti¬ 
cal studies show that the intensity of the produced mag¬ 
netic field rises approximately linearly with the centre of 
mass energy (-^/snn) of the colliding nucleons [1, 2]. The 
Lorentz boosted electric fields in such collisions also be¬ 
comes very strong which is same order of magnitude as 
magnetic field {eB ^ eE ^ lOm^ for a typical An-An col¬ 
lision at top RHIC energy -\/5 nn= 200 GeV), where rrij^ is 
the pion mass. Such intense electric and magnetic fields 
are strong enough to initiate the particle production from 
vacuum via Schwinger mechanism [3] . Using quantum 
chromodynamics it was shown in Ref. [4] that beyond a 
critical value of magnetic field the quark-antiquark state 
can possibly attain negative mass (in the limit of large 
number of colours). Thus it is important to know if there 
is truly an upper limit of magnetic field intensity allowed 
by the quantum chromodynamics when applied in heavy 
ion collisions. Or the magnetic field can grow to arbi¬ 
trary large value with increasing y^SNN as predicted in 
some earlier studies [1,2]. In this work we shall calculate 
the electromagnetic field intensity without considering 
any such constraints, i.e, we assume that the electric and 
magnetic fields can attain any arbitrary large values. 

There are several other interesting recent studies re¬ 
lated to the effect of ultra-intense magnetic fields in 
heavy-ion collisions. Here we briefly mention a few of 
them which might be relevant to the present study. In 
presence of a strong magnetic field as created in heavy- 
ion collisions, a charge current is induced in the Quark 
Gluon Plasma (QGP), leading to what is known as the 
“chiral magnetic effect” (CME) [5]. Within a 3+1 di¬ 
mensional anomalous hydrodynamics model a charge de¬ 
pendent hadron azimuthal correlations was found to be 
sensitive to the CME in Ref. [6]. Along with CME, it 
was also predicted theoretically that massless fermions 


with the same charge but different chirality will be sep¬ 
arated, yielding what is called the “chiral separation ef¬ 
fect” (CSE). A connection between these effects and the 
Berry phase in condensed matter was also pointed out in 
Ref. [7-9]. In hadronic phase a significant changes in the 
hadron multiplicity was observed in presence of a strong 
magnetic field within a statistical hadron resonance gas 
model in Ref. [20]. There are lot of other important 
relevant work in this new emerging field which we can¬ 
not refer here, one can see recent reviews on this topic in 
Ref. [10-12] for more details. 

The relativistic hydrodynamic models have so far 
nicely explained the experimentally measured anisotropic 
particle production in the azimuthal directions in heavy 
ion collisions. The success of hydrodynamics model 
shows that a locally equilibrated QGP with small ratio 
of shear viscosity to entropy density is formed after the 
collision within a short time interval ^0.2-0.6 fm [13-19]. 
It is also well known that the final momentum anisotropy 
in hydrodynamic evolution is very sensitive to the initial 
(geometry) state of the nuclear collisions. So far almost 
all the hydrodynamic models studies have neglected any 
influence of magnetic fields on the initial fluid energy- 
density or on the space-time evolution of QGP. But as 
we know the initial magnetic field is quite large, it is 
important to investigate the relative importance of large 
electro-magnetic field on the usual hydrodynamical evo¬ 
lution of QGP. Eor that one need a full 3+1 dimensional 
magnetohydrodynamic code to numerically simulate the 
space time evolution of QGP with magnetic fields. While 
one can gain some insight about the relative importance 
of the magnetic filed on the initial energy density of 
the QGP fluid by estimating the quantity plasma sigma, 
which is the dimensionless ratio of magnetic field energy 
^ to the fluid energy density( 5 ) : cr = ^. In plasma 
a indicates that one can no longer neglect the ef¬ 
fect of magnetic fields in the plasma evolution (in some 
situation a ~ 0.01 may also effect the hydrodynamic evo¬ 
lution) [21-24]. In the present study we use MC-Glauber 
model [25, 26] to calculate e-by-e magnetic fields and fluid 
energy density in An-An collisions at = 200 GeV 
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and investigate the relative importance of the magnetic 
field on initial fluid energy density. 

As mentioned earlier, the typical magnetic field pro¬ 
duced in a mid-central Au-Au collisions at = 

200 GeV reaches ^ lOm^, which corresponds to field en¬ 
ergy density of ^ 5GeV/fm^. Hydrodynamical model 
studies show that the initial energy density for such cases 
is ^ lOGeV/fm^, thus implying a ^ 0.2 under these 
conditions. However, the magnetic field produced at the 
time of collisions decays very quickly if QGP does not 
possess finite electrical conductivity [27-29]. Thus in or¬ 
der to correctly estimate cr, one need to consider the 
proper temporal evolution of magnetic fields until the 
thermalisation time (r^ ^ 0.5/m for Au-Au collisions at 
RHIC) when the hydrodynamic evolution starts. Since 
the spatial distribution of fluid energy density as well as 
the electromagnetic fields varies e-by-e we also calculate a 
accordingly. The spatial distribution of electric and mag¬ 
netic fields in heavy ion collisions was previously studied 
in Ref. [30, 31]. 

In present work we study the spatial distribution of a 
in y/^NN = 200 GeV Au-Au collisions for two different 
impact parameters (b=0, and 12 fm). The temporal evo¬ 
lution of the magnetic fields after the collision is taken 
into account in a simplified manner which will be dis¬ 
cussed in the next section. We also investigate the cor¬ 
relation between a and fluid energy density in the trans¬ 
verse plane. The paper is organised as follows: in the 
next section, we discuss about the formalism. Our main 
result and discussion are presented in section III. A sum¬ 
mary is given at the end in section IV. 


II. FORMALISM AND SETUP 


We constructed a spatial grid of size 10 fm in each di¬ 
rection (X and Y) with the corresponding grid spacing of 
Ax = Ay = 0.5 fm for e-by-e calculation of electromag¬ 
netic fields and fluid energy density in the plane trans¬ 
verse to the trajectory of the colliding nuclei. The posi¬ 
tion of colliding nucleons are obtained from MG-Glauber 
model in e-by-e basis. The position of nucleons are ran¬ 
domly distributed according to the Wood-Saxon nuclear 
density distribution (as shown in Fig ( 1 )). We adopt 
the usual convention used in heavy ion collisions for de¬ 
scribing the geometry of the nuclear collisions , i.e., the 
impact parameter vector ( 6 ) of the collision is along X 
axes and the colliding nuclei are symmetrically situated 
around the (0,0) point of the computational grid. The 
electric and magnetic fields at point r(x, y) at time t due 
to all charged protons inside two colliding nucleus are 
calculated from the Lienard-Weichart formula 


Nr, 




2=1 


Ri Ri Vi 

^Ri Ri ' ^2^ 
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Nproton 
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Where E and B are the electric and magnetic field vec¬ 
tor respectively, e is the charge of an electron, Z is the 
number of proton inside each nucleus, Ri = x — Xi{t) is 
the distance from a charged proton at position Xi to x 
where the field is evaluated, Vi is the velocity of the Tth 
proton inside the colliding nucleus. Ri is the magnitude 
of Ri. The summation runs over all proton {N proton) 
inside the two colliding nuclei. Following Ref. [I] we as¬ 
sume that because of the large Lorentz factor (7 ^ 100) 
the colliding nuclei are highly Lorentz contracted along 2 : 
direction and all the colliding protons have same velocity 

^first _ pjSecond _ 

lated to the c.m energy (-^/^nn) through the relationship 

Vz = ~ (’ where rup is the proton mass. Note 

that according to Eq. (I) and (2) the electric and mag¬ 
netic fields diverge as ^ 0 , to remove this singularity 
we assume a lower value Rout = 0.3 fm as used in Ref [I]. 
This particular value of Rout = 0-3 fm was fixed as an 
effective distance between partons and it was found that 
the calculated electromagnetic field has weak dependence 
for 0.3/m < Rout ^ 0.6/m. We note that the quantities 
eB and eE has dimension GeV^ and the conversion from 
GeV^ to Gauss is given by IGeV^ = 5.12 x 10 ^^ Gauss. 
It is customary to use Milne co-ordinate (r = 

y /{E — 2 :^),X, T] = in heavy ion collisions. 

For our case we shall concentrate on the mid-rapidity 
region (z^O) where t ^ r. 

By using the MG-Glauber model we also compute the 
fluid energy density in transverse plane from the po¬ 
sition of wounded nucleons. This is a common prac¬ 
tice to initialise energy density for e-by-e hydrodynamics 
simulations. Since the positions of the wounded nucle¬ 
ons (Vwound) are like delta function in co-ordinate space, 
in order to calculate the energy density profile for hy¬ 
drodynamics simulations one need to smooth the initial 
profile by introducing Gaussian smearing for every collid¬ 
ing nucleons. The fluid energy density 5 is parameterised 
as 


/ N ^ou nd (^x-Xi(b))‘^ + {y-yi (b))^ 

e (^x,y,ag,bj = k ^ e , (3) 

here x, y is the co-ordinate of computational grid, 
Xi(b),yi{b) are the co-ordinate of wounded nucleons for an 
impact parameter 6 , ag is the Gaussian smearing which 
is taken to be 0.5 fm (unless stated otherwise) for our 
calculation, /c is a constant which is tuned to match the 
initial central energy density for event averaged case. We 
estimate k=6 which results in the initial central energy 
density 40GeV//m^ for 5 = 0 fm collision. This is the 
typical value of initial energy density used in e-by-e hy¬ 
drodynamics model to reproduce the experimental mea- 
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sured charged particle multiplicity in Au-Au collisions 
at ysNN =200 GeV for an initial time ^ 0.5 fm [16]. 
The same k factor is used to calculate the initial energy 
density for all other impact parameter. 

Once we calculate the electromagnetic field and the 
fluid energy density in the transverse plane, the plasma 

a{x,y,b) = ^ readily obtained for each event. 

For our case we only considered the transverse compo¬ 
nents Bx and By to calculate the total magnetic energy 
density, since B^ B^^By. As mentioned, the hydro¬ 
dynamics expansion of the QGP fluid starts after a time 
^0.5 fm, and because of the relativistic velocities of the 
charged protons the produced magnetic fields decays very 
quickly. If there is no conducting medium then the mag¬ 
netic field decay as ^ . But in presence of a conduct¬ 

ing medium the decay can be substantially delayed [29]. 
However, the thermodynamic and transport properties 
of the nuclear matter right after the collision upto the 
time when the system reaches local thermal equilibrium 
is poorly known. Thus we investigate in our study two 
different scenarios when calculating a{x,y,b). From now 
on we will omit b in the expression of a, and because of 
spherical symmetry of the colliding nuclei we omit the 
vector arrow and simply write b for the impact parame¬ 
ter. 

(i) In the first scenario, following Ref. [29] we assume 
that the matter in pre-equilibrium phase has finite elec¬ 
trical conductivity and the field components are evalu¬ 
ated at Ti =0.5 fm by reducing the magnitudes of the 
initial magnetic field (at r=0 fm) to O.I times. This is a 
simplification of the actual scenario, since the time evo¬ 
lution of the fields depend on the electrical conductivity, 
the impact parameter and on the Lorentz gamma( 7 ) of 
the collisions. According to Ref. [29] the initial electro¬ 
magnetic field produced in a b=7 fm collision and for an 
electrical conductivity crez=5.8 MeV reduce ^50% to its 
original value after r ^0.5 fm. Note that for simplifi¬ 
cation in the numerical simulation we have ignored the 
impact parameter dependence of the evolution of electro¬ 
magnetic field in medium as was discussed in Ref. [29]. 

(ii) In the second scenario, we assume the magnetic 
field is evolved in vacuum (zero electrical conductivity) 
until the hydrodynamics expansion starts. For this case 
we reduced the magnitude of the initial electromagnetic 
field 0.01 times. 

We note that in reality the situation may he in be¬ 
tween the above mentioned two scenarios. From now on 
we denote the first and second scenario by medium and 
vacuum respectively. 

We consider 1000 nucleus-nucleus collisions for our 
present calculation for each impact parameter. 


III. RESULTS AND DISCUSSION 

At first we shall concentrate on the electromagnetic 
fields computed at the centre of the fireball (i.e. at point 
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FIG. I: Distribution of nucleons inside target and projectile 
nuclei in a typical Au-Au collisions at ^snn =200 GeV for 
b=I2 fm. 
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FIG. 2: Impact parameter dependence of event averaged mag¬ 
netic and electric fields at the centre of the fireball for Au-Au 
collisions at Vsnn= 200 GeV. 


X = y = 0 in our computational grid). Fig. 2 shows the 
event averaged value of magnetic and electric fields as a 
function of impact parameter b. The By , its absolute 
value \By\, and x component of the electric field Ex are 
shown by pink dashed, blue dotted, and black solid lines 
respectively. We note that our result is consistent with 
the result of Ref. [I]. We also checked other components 
of electric and magnetic fields and they are found to be 
consistent with Ref. [I]. 

The electric and magnetic fields are created in high 
energy heavy-ion collisions in presence of the electri¬ 
cally charged protons inside the two colliding nucleus. 
Whereas both neutron and protons inside the colliding 
nuclei deposit energy in the collision zone as a result of 
elastic and inelastic collisions among them. Since the po¬ 
sitions of protons in the colliding nucleus is different with 
that of the positions of all nucleons, resulting spatial dis¬ 
tribution of the electromagnetic field is expected to differ 
from that of the initial fluid energy density. Fig. 3 shows 
the event averaged value of fluid energy density at point 
(x = y = 0) as a function of impact parameter b. The en- 
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FIG. 3: Impact parameter dependence of event averaged cen¬ 
tral energy density (< £ (0, 0) >) of fluid for An-An collisions 
at =200 GeV. 


ergy density is obtained from Eq. 3 for k =6. This specific 
value of k was chosen in order to obtain the central energy 
density ^ 40 GeV/fm^ for 6 = 0 fm collisions. From pre¬ 
vious studies [16] we note that the initial central energy 
density for central (0 — 5% centrality which corresponds 
to 6 ^ 2 fm) An-An collisions requires 5 ^ 40 GeV/fm^ 
at initial time = 0.6 fm at the centre of the fireball 
(x = y = 0) to reproduce the experimentally measured 
charged hadron multiplicity at y^^NN = 200GeV. How¬ 
ever, we note that a different initial time (r^) will give 
different initial energy density [32] , in that case the mag¬ 
nitude of magnetic field at will also be different. From 
Fig. 2 and Fig. 3 we notice that fluid energy density de¬ 
creases whereas the intensity of magnetic field increases 
with b. It is thus expected that a{x,y^b) will reach its 
maximum value for b ^ 12 fm. So far we have shown the 
event average E{x^y)^B{x^y) and e{x,y) at the center 
of the collision zone. 

Top panel of Fig. 4 shows the event averaged e{x^y) for 
b =0 fm collisions. Since the An nucleus is almost spher¬ 
ical in shape, a head on An-An collision deposits energy 
in an almost circular zone. Different colour schemes in 
the legend denotes the energy density in unit of GeV/fm^. 
Middle and bottom panel of Fig. 4 shows the correspond- 

73 2 

ing magnetic field energy density due to the y and x 
component of B respectively, where B is calculated at 
time r = 0. We observe that the distribution of mag¬ 
netic field energy is similar to the fluid energy density 
obtained from elastic and inelastic nucleon-nucleon colli¬ 
sions in MC-Glauber model. The magnetic field energy 
density due to Bx and By is peaked at the centre and has 
a SO(2) rotational symmetry for b=0 fm collision. This 
is not surprising since the positions of the protons for 
b =0 fm collisions have such rotational symmetry about 
the centre of the fireball in the transverse plane for event 
averaged case. The situation for a non-zero impact pa¬ 
rameter collisions becomes different . The overlap zone 
between the two nuclei becomes elliptical, as can be seen 
from the top panel of Fig. 5 which corresponds to e{x^ y) 
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FIG. 4: Top Panel: 1000 event averaged initial energy den¬ 
sity of QGP for b=12 fm An-An collisions at V'®nn= 200 GeV. 
Middle Panel: 1000 event averaged magnetic field energy den¬ 
sity calculated from y component of the magnetic field for 
b=12 fm An-An collisions at y^NN GeV. Bottom Panel: Same 
as middle panel but for the x component of the magnetic fields 
Bx- 


for b =12 fm. The middle and bottom panel of Fig. 5 
shows the corresponding energy density for By and Bx 
components. We find that the field energy density due 
to By has similar shape as fluid energy density, but that 
due to Bx has maximum in a dumbbell shaped region 
which is different from the initial fluid energy density. 

So far we have shown event averaged value of 5 and 
components of H. It is not clear from the above discus¬ 
sion whether the magnetic field energy density is negligi¬ 
ble compared to the initial fluid energy density for every 
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FIG. 5: Top Panel: 1000 event averaged initial energy density 
of QGP for An-An at y^SNN =200 GeV for impact parameter 
b=12 fm collisions. Middle Panel: event averaged magnetic 
field energy density calculated from y component of the mag¬ 
netic field for b=12 fm An-An collisions at y^SNN GeV. Bot¬ 
tom Panel: Same as middle panel but for the x component of 
the magnetic fields Bx. 
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events because both 5 and are lumpy in the transverse 
plane as shown in Fig. 6. This leads us to study cr{x,y) 
on e-by-e basis. 


A. Event-by-event a(x,y) 

Top panel of Fig. 6 shows the energy density, middle 
and bottom panels show corresponding a (x, y) at r =0.5 
fm for evolution of the magnetic field in medium and in 
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FIG. 6: Top panel: fluid energy density , Middle panel: 
a {x,y) for the medium , Bottom panel: a {x,y) for the vac¬ 
uum in a single b=12 fm An-An ^snn= 200 GeV collision. 
The shaded annular region in middle and bottom panel cor¬ 
responds to 0.01 < a {x,y) < 10 . 


vacuum respectively for a single event of b=12 fm colli¬ 
sions. The shaded band in the middle and bottom panels 
correspond to the zones where 0.01 < a{x,y) < 10 (in¬ 
creasing in the outward direction). As expected a {x,y) 
reaches its maximum value in regions where e{x^y) be¬ 
comes small. However, note that those regions of large 
a {x,y) strongly depends on the temporal evolution of 
the magnetic field from r = Ofm until the hydrodynam¬ 
ics expansion starts at time r^. This can be seen from 
the bottom panel of the same figure where the regions 
of large a {x,y) moves outward as the magnetic filed for 
this case decays faster than the case of medium with fi¬ 
nite electrical conductivity. 

We observe here that even if the magnetic field de¬ 
cays quickly (as in vacuum) until the hydrodynamics 
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FIG. 7: Top panel: e-by-e distribution of cr (0, 0) as a function 
of £ (0,0) for b=0 fm An-An collisions at ^snn =200 GeV. 
Bottom panel: same as top panel but for (x = 7, y = 0). 


expansion starts, there is a corona of large cr(x,^). It 
is then important to consider magnetohydrodynamics 
framework to investigate further the possible effects of 
those large a (x, y) zone on the space-time evolution of 
the QGP ffuid.We expect that since the region of large 
cr seems to lie mostly outside the places where e{x^y) is 
high there will be minor modification in the transverse 
evolution of the QGP fiuid when the effect of magnetic 
field is taken into account. The above conclusion is made 
by investigating only one particular event, in order to un¬ 
derstand the ensemble of events let us look at the e-by-e 
distribution of cr (x, y) at the centre (x = y = 0) and at 
the periphery of the collision zone. 

Top panel of Fig. 7 shows the e-by-e distribution of 
cr (0,0) as a function of 5 (0, 0) for b=0 fm collisions. The 
bottom panel of the same figure shows the event distribu¬ 
tion of cr (7, 0) versus 5 (7, 0). All results are obtained for 
magnetic field evolution in medium. Naively one expects 
that a{x,y) and e{x,y) should be anti-correlated, i.e., 
for places where 5 is large cr will be small and vice-versa, 
the same conclusion was made in Ref. [31]. But it is clear 
from Fig. 7 that there is no such simple relationship be¬ 
tween 5 and cr for b=0 fm collisions in MC-Glauber model 
. In fact, we notice that at the center of the collision zone 
5 and cr are almost uncorrelated. For regions at the pe¬ 
riphery of the collision zone (bottom panel) we observe 
similar behaviour, but notice that here cr (7, 0) may reach 
^ 1 in some events, whereas for x = ^ = 0 it never exceed 
0 . 01 . 

Now let us discuss the result for b=12 fm collisions 
where the relative importance of magnetic field is ex- 
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FIG. 8: Top panel: e-by-e distribution of a (0, 0) as a function 
of £ (0,0) for b=12 fm An-An collisions at ^snn =200 GeV. 
Bottom panel: same as top panel but for (x = 3, y = 0). 

pected to be highest. Top panel of Fig. 8 shows the e-by- 
e distribution of cr (0, 0) as a function of 5 (0, 0) for b=12 
fm collisions. We notice that like cr (0, 0) distribution for 
b=0 fm collisions, most of the events have cr (0, 0) ^ 0.01. 
However, for few events cr(0,0) ^ 1. Like b = Ofm col¬ 
lisions here we also notice no clear correlation between e 
and cr. The bottom panel of Fig. (8) shows the distribu¬ 
tion of cr (3,0) as a function of 5 (3, 0). We notice that a 
considerable number of events have cr ^ 1 for this case. 

Next we discuss the event averaged transverse profile 
of a{x^y) for b=0 and 12 fm as depicted in Fig. 9 and 
10 respectively. As expected, the event averaged a {x,y) 
in the range 0.01 < cr < 10 for b=0 fm collisions (Fig. 9) 
form an annular region enclosing the high energy density 
zone of the QGP fiuid. Top panel of Fig. (9) shows the 
result for magnetic field evolution in vacuum and the 
bottom panel shows in medium. However, cr (x, y) for 
b=12 fm collisions shows different spatial distribution as 
depicted in Fig. 10. The non-trivial contour in this case 
results from the fact that for some events cr (x, y) becomes 
very large and the event averaged value is dominated by 
those large cr. The top and bottom panel shows the result 
for vacuum and medium respectively. 


B. Sensitivity of a(x,y) on Gaussian smearing 

The Gaussian smearing cr^ in Eq. (3) is a free param¬ 
eter which is usually taken in the range ^ 0.1-1.0 fm. 
Here we discuss the sensitivity of our result on Gaussian 
smearing by setting cr^ = 0.25fm which is taken from a 
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FIG. 9: Event averaged a {x, y) in the range 0.01 < a < 10 
(shaded region) for Au-Au collisions of b=0 fm at = 200 

GeV . Top panel: for vacuum. Bottom panel: for medium. 


10 “ 

10-1 
10-2 
D 10-3 
10-4 
10 -“ 

10-6 
10-2 
10 “ 

10-1 
10-2 
° 1 (H 
10-4 
10-5 
10-6 
10-2 

10-5 10-2 10-1 10 “ 10 ' 102 

e(GeV/fm®) 

FIG. 11: Top panel: e-by-e distribution of cr (0, 0) as a func¬ 
tion of £ (0, 0) for Au-Au b=0 fm collisions at ^snn =200 
GeV. Bottom panel: same as top panel but for {x = 7,y = 0). 
(7j= 0.25 for both the cases. 





FIG. 10: Event averaged a {x, y) in the range 0.01 < f < 10 
(shaded region) for Au-Au collisions of b=12 fm at ^snn = 
200 GeV. Top panel: for vacuum. Bottom panel: for medium. 
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FIG. 12: Top panel: e-by-e distribution of a (0, 0) as a func¬ 
tion of £(0,0) for Au-Au b=12 fm collisions at y^SNN =200 
GeV. Bottom panel: same as top panel but for (x = 3, y = 0). 
crq=0.25 for both the cases. 
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recent study [33]. Reducing (jg results in much lumpy 
initial energy density hence we expect a different spa¬ 
tial dependence of cr(x, y) compared to the previous case 
where (jg = 0.5fm is used. For (jg = 0.25fm we adjusted 
k to 8i new value k=17 to keep the event-averaged ini¬ 
tial central energy density for b=0 fm collisions same as 
before i.e., ^ 40GeV/fm^ . Top panel of Fig. 11 shows 
the e-by-e distribution of cr (0,0) as a function of 5 (0, 0) 
for b=0 fm An-An collisions. Bottom panel shows the 
same but for x = 7,y = 0. Comparing Fig. 7 and 11 we 
found that changing ag from 0.5 fm to 0.25 fm changes 
the e-by-e distribution of a vs 5. Since the energy den¬ 
sity is more lumpy for (jg = 0.25 fm than 0.5 fm, the 
number of events with large a increases. To see the ef¬ 
fect of changed (jg in peripheral collisions, we show the 
e-by-e distribution of a vs 5 for b=12 fm in Fig. 12. Top 
panel of Fig. 12 shows the e-by-e distribution of a (0, 0) 
vs 5 (0,0) and the bottom panel shows e-by-e distribu¬ 
tion of cr(3,0) vs 5(3,0). It is clear that for b=12 fm 
collisions the correlation between a and 5 at the center 
(x = ^ = 0) is sensitive to cr^, and the maximum value 
of cr ^ 1, in contrary to what was observed for the case 
b=0 fm collisions. 


IV. SUMMARY 

We have studied the relative importance of magnetic 
field energy on initial fiuid energy density of the QGP by 
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evaluating a = ^ for An-An collisions at ■s/s^= 200 
GeV. The fiuid energy density and electromagnetic fields 
are computed by using MC-Glauber model. The electro¬ 
magnetic field and initial fiuid energy density are calcu¬ 
lated by using following parameters: the cutoff distance 
Rent = 0.3fm, Gaussian smearing parameter cTg =0.5, 
(and 0.25 fm) and the scalar multiplicative factor k =6, 
(and 17). The initial energy density (at time = 0.5 
fm) for the fiuid is fixed to ^ 40 GeV/fm^. The ratio 
of the magnetic field energy density to the fiuid energy 
density a is evaluated in the transverse plane for two dif¬ 


ferent impact parameters b =0 , and 12 fm. We find 
that for most of the events, at the centre of the colli¬ 
sion zone a (0, 0) ^ 1 for both b=0, and 12 fm collisions. 
However, at the periphery of the collision zone where 
5 becomes small we observed a region of large a. For 
large impact parameter collisions a becomes larger for 
peripheral collisions (large b) compared to central (small 
b) collisions as a result of increase in magnetic field and 
decrease in fiuid energy density. We observe that in cen¬ 
tral collisions (b= 0 fm) at the center of collision zone 
a for most of the events. However, large a is ob¬ 
served in the outer regions of collision zone. In periph¬ 
eral collisions a becomes quite large at both center and 
periphery of the collision zone. From this observation 
we conclude that initial strong magnetic field might con¬ 
tribute to the total initial energy density of the An-An 
collisions (or other similar heavy ion collisions like Pb- 
Pb) significantly. However, the true effect of large a (or 
large magnetic fields) will remain unclear unless one per¬ 
forms realistic magneto-hydro dynamics simulation with 
the proper initial conditions, for example see Ref. [21- 
24] for some theoretical estimates. Note that the result 
in this paper are obtained for a specific model of initial 
conditions (MC-Glauber model) with few free parame¬ 
ters. We have not explored all possible allowed values 
of these free parameters. In future we may incorporate 
other initial conditions and more realistic time evolution 
of the electromagnetic fields in the pre-equilibrium phase 
(as described in Ref. [29]) to study the effect of mag¬ 
netic fields on initial fiuid energy density distribution. It 
is also interesting to study the similar thing for lower 
collisions where the decay of magnetic field in vac¬ 
uum is supposed to be much slower than the present case 
because of the slower speed of colliding nuclei, and also 
the corresponding initial energy density for such cases 
is smaller than the present case of An-An collisions at 
-^/snn = 200GeV. 
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